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Abstract 



>. 

CN . We compute the tunneling probability from the symmetric phase to the true 

t — ■ 

£NJ ■ vacuum, in the first order electroweak phase transition of the MSSM, and the 

q ' corresponding Higgs profiles along the bubble wall. We use the resummed 

two-loop temperature-dependent effective potential, and pay particular at- 
^ e— I tention to the light stop scenario, where the phase transition can be suffi- 

ciently strongly first order not to wipe off any previously generated baryon 



asymmetry. We compute the bubble parameters which are relevant for the 
baryogenesis mechanism: the wall thickness and A/3. The two-loop correc- 
tions provide important enhancement effects, with respect to the one-loop 
results, in the amount of baryon asymmetry. 
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1 Introduction 



Electroweak baryogenesis |1] is an appealing mechanism to explain the observed value 
of the baryon-to-entropy ratio, tib/s ~ 10~ 10 , at the electroweak phase transition [0, [|, 
that can be tested at present and future high-energy colliders. Although the Standard 
Model contains all the necessary ingredients JjJ for a successful baryogenesis, it fails 
in providing enough baryon asymmetry. In particular it has been proven by perturba- 
tive ||, |, U and non-perturbative |7j methods that, for Higgs masses allowed by present 
LEP bounds, the phase transition is too weakly first order, and any previously gener- 
ated baryon asymmetry would be washed out after the phase transition. On the other 
hand the amount of CP violation arising from the CKM phase is too small for gener- 
ating the observed baryon asymmetry ||. Therefore electroweak baryogenesis requires 
physics beyond the Standard Model at the weak scale. 

Among the possible extensions of the Standard Model at the weak scale, its minimal 
supersymmetric extension (MSSM) is the best motivated one. It provides a technical 
solution to the hierarchy problem and has its roots in more fundamental theories unify- 
ing gravity with the rest of interactions. The MSSM has new violating phases || that 
can drive enough amount of baryon asymmetry [[H]]-[|nj provided that the previous 
phases are not much less than 1 and the charginos and neutralinos are not heavier than 
200 GeV []. As for the strength of the phase transition [17|]-[19]' a re gi° n i n the space 
of supersymmetric parameters has been found [^0J-p8[ where the phase transition is 
strong enough to let sphaleron interactions go out of equilibrium after the phase tran- 
sition and not erase the generated baryon asymmetry. This region (the so-called light 
stop scenario) provides values of the lightest Higgs and stop eigenstates which are: 75 
GeV < rrih < 105 GeV, 100 GeV< m t . It will be covered at LEP2 and Tevatron 

colliders. 

In all calculations of the baryon asymmetry the details of the wall parameters play a 
prominent role in the final result. In particular the wall thickness, L u , and the relative 
variation of the two Higgs fields along the wall, A/3, are typical parameters which the 
generated baryon asymmetry depends upon. Although reasonable assumptions about 
the Higgs profiles along the wall have been done, as e.g. kinks or sinusoidal patterns 
interpolating between the broken and the symmetric phases, as well as estimates on the 
value of A/3 based on purely potential energy considerations, it is clear that the relia- 
bility of those estimates as well as more precise computations of the baryon asymmetry 



1 The first and second generation squarks are required to have masses O(few) TeV because of the 
experimental bounds on the neutron electric dilope moment. 



1 



should rely on realistic calculations of the Higgs profiles and the tunneling processes 
from the false to the true vacuum. Such a task, having been achieved in the case of one 



Higgs field in the Standard Model ||29|| , was still missing in the case of two-Higgs field 
models. 

In this paper we compute the tunneling probability from the false to the true vacuum 
in the first order electroweak phase transition of the MSSM, and the corresponding 
Higgs profiles along the wall. In particular we will concentrate in the region of the 
supersymmetric parameters corresponding to the light stop scenario, where the phase 
transition is strong enough not to wash out the generated baryon asymmetry. We will 
use the MSSM effective potential at finite temperature including the most important 
two-loop corrections. The two-loop corrections have been proven to be very important 
in determining the strength of the phase transition and we will demonstrate that they 
also play a very relevant role in the baryogenesis mechanism. In particular we will 
compare the results using the one-loop effective potential, that has been used so far for 
the determination of the baryon asymmetry, with the corresponding two-loop results. 
We will see that there is an important enhancement coming from two-loop effects. 
The plan of the paper is as follows. In section 2 the method to compute the bubbles 
in the MSSM will be described in detail. The results of our numerical analysis will 
be presented in section 3. The possibility of bubbles involving squark fields will be 
explored in section 4 and section 5 is devoted to draw our conclusions. 



2 Higgs Bubbles 



The Higgs sector of the MSSM requires two Higgs doublets, with opposite hypercharges, 
as 

/ \ / 




-1/2 \ 2 / 1/2 



and tree-level potential: 



V {0) = mlHlHt + mlHlH 2 + m^(#i# 2 + /i.e.) (2.2) 



The field configuration describing the tunneling in a theory with just one scalar field 
in three or more dimensions is known to be spherically symmetric [ 3"0fl . Here we will 
assume that the bubbles driving the electroweak phase transition in the MSSM and 
involving the two neutral higsses also have spherical symmetry. They correspond to the 



2 



ansatz: 



Hi{x) = hi{r) 



1 





H 2 {x) = h 2 {r) 




1 



(2-3) 



where r = yx 2 . In the presence of the background (p.3|) the tree-level potential reads 



as 



V^[h u h 2 ) = m\ h\{r) + m\ h 2 2 (r) + 2 m^i(r)/i 2 (r) + 



9 z + 9 



12 



h\{r)-h 2 2 {r) (2.4) 



and the one-loop corrections, in the 't Hooft-Landau gauge and in the MS renormaliza- 
tion scheme, are given by 



V^(h 1 ,h 2 )=Y / -^- 2 mt(h 1 ,h 



log 



m 2 (h 1: h 2 



a 



(2.5) 



where Q is the MS renormalization scale, Cj=3/2 (5/6) for scalar bosons and fermions 
(for gauge bosons), m 2 (hi,h 2 ) is the field dependent mass of the i th particle in the 
background hi, h 2 , and n« is the corresponding number of degrees of freedom, which is 
taken negative for fermions. 

By minimizing the effective potential + with respect to (hi,h 2 ), and im- 
posing the minimum of the potential to be at (v i, v 2 ), with v = yv \ + v \ = 174.1 GeV, 
and tan/5 = v 2 /vi fixed, we can eliminate m\ and m 2 in favour of the other parameters 
of the theory, while m| can be traded in favour of the one-loop corrected squared mass 



m A of the CP-odd neutral Higgs boson ||19|| . 



Then we can write the finite temperature effective potential as 

V(hi, h 2 , T) = V^(hi, h 2 ) + VW{h u h 2 ) + AV(hi, h 2 , T) 



(2.6) 



where the thermal correction AV contains the one- loop pf, AV^{h x ,h 2 ,T), plus 
the leading two-loop radiative corrections [^5], ^6], p8| , AV^ 2 >{h\, h 2 , T), of the daisy 
resummed theory. Since the two-loop corrections have been proved to be very relevant 
for the strength of the phase transition |25j , and present estimates of bubble parameters 
are based on one-loop thermal corrections, we will often compare our numerical two-loop 
calculations of these parameters with one-loop results. 
The Euclidean action of configuration fl2.3|) is given by: 



S 3 (T) =4njdr r 2 



' d h\ 



d ho 



V{hi,h 2 ,T) 



(2.7) 



dr J \dr J 

where V(hi, h 2 , T) is the effective potential given in fl2.6Q , shifted in such a way that 



3 



V(0, 0, T) = 0. The bubble is then the solution of the equations 

d 2 h x 2 dh x 1 dV 
+ -- 



dr 2 r dr 2 dh 

supplied by the boundary conditions 



dr 2 r dr 2 dhi 

d 2 h 2 2 dh 2 1 dV 



dhi 
dr 



= 0; ^Uo = (< = 1,2). (2.9) 

r=0 

This solution only exists in the range T d < T < T , where T d is the temperature at 
which the minimum at (hi, h 2 ) is degenerate with that at the origin, and T is the 
temperature at which the origin gets destabilized along some direction and becomes a 
saddle point. 

Solving the previous differential equations is a difficult task. There is no systematic 
procedure to find the solution. For example, the usual overshooting-undershooting 
method cannot be implemented in the two (or more) field case. A possible technique is 
to define a functional, acting in some space, having its minimum at the solution of the 
previous equations. Notice that the bubble solution is a saddle point of the Euclidean 
action, not a minimum, so we cannot use directly S3. One possibility, explored in 
reference [ftl], is to add to S3 some (nonlocal) terms that lift the falling directions. 



Then the bubble is a minimum of this improved action. Here we will follow a slightly 
different approach. We will consider the following functional f| 



T = 47r dr r 



' d 2 h x 2 dh x 
+ 



1 dV 



dr 2 



r dr 2 dh\ 



+ 



'fh 
dr 



2 dho 

+ 



1 dV 



r dr 2 dhi 



(2.10) 



By construction, the bubble solution is a minimum of the functional T. In order 
to find numerically this solution, we have discretized the radial integral and spatial 
derivatives in the functional. The length scale used for the radial variable is M^} 1 . 
The bubble typically spreads over a range of order 100 M^- Instead of minimizing 
simultaneously 2N points describing the hi and h 2 profiles on the lattice, we will use 
an iteration method. First, to have an estimate of the bubble size, we solve a reduced 
problem. Assuming that the tunneling takes place along the direction given by con- 
stant tan/3(T), we can use the standard overshooting-undershooting method to find 
numerically this configuration. This provides us with a good estimate of the size of the 



2 This kind of functionals has also been used to obtain unstable solutions, such as sphalerons, on 



the lattice 32 
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lattice, where to solve the discretized problem, and also a first starting configuration 
for the minimization. Now, we fix one of the Higgses, hi for example, to the previous 
estimate and allow the N points describing the other Higgs to change, minimizing the 
functional given in (|2.10|) . We iterate this procedure, now fixing the h 2 field. We have 
used iV = 140 for the discretization. In the range of the supersymmetric parameters 
we have explored, we found a good convergence after two or three iterations in hi, h 2 . 
Various tests derived from the fact that the obtained solutions are extremals of the 
Euclidean action, as e.g. the virial theorem, have been proved numerically to hold. We 
have also verified stability of the profiles against increasing values of N. 



3 Numerical results 

The strength of the phase transition, after introduction of the leading two-loop thermal 



corrections, has been studied in Refs. [25]- [28]. The light stop scenario was carefully 



analyzed in Ref. [£8|] where the effective potential along both the Higgs and £r di- 
rections was considered, including the corresponding two-loop corrections. There, the 
different possibilities according to the cosmological evolution of the fields were classi- 
fied as regions of stability, metastability and instability of the Higgs minimum, and 
regions where a two-step phase transition can proceed. In this paper we concentrate 
specifically in the stability region, since the metastability and two-step regions would 



require considering the tunneling to (and from) the color breaking minimum [33], which 



is outside the scope of the present paper. We will closely follow the allowed parameter 



space found in Ref. [2S 



We choose tuq = 1 TeV in such a way that the supersymmetric corrections to 
the p-parameter become small, giving hence a good fit to the electroweak precision 
data coming from LEP and SLC On the other hand we will take as reference point 
At = At — \ij tan/3 = 0, mj = 150 GeV and tan/3 = 2.5, itla = 200 GeV, providing 
rrih = 80 GeV, which is consistent with present experimental bounds on the MSSM 
lightest Higgs mass. This particular point in the plane (mh,Tn~) is acceptable from 
the point of view of avoiding the wipe off of the baryon asymmetry after the phase 
transition, as was shown in Figs. 1 and 2 of Ref. [p8"| . We have checked that the 
obtained results are rather generic, as we will explicitly show by varying the parameter 
rriA, so for other points in the allowed region of the plane (m^, m-r) we obtain similar 
features. 

In Fig. 1 we plot the euclidean action S3 ( |2.7| ) for the bounce solution (|2.8|) as a 



function of the temperature for A t = 0, m j = 150 GeV, tan/3 = 2.5 and tua = 200 
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GeV. We can see from Fig. 1 that S3 goes to infinity at the degeneracy temperature 
Td ~ 96 GeV, while it goes to zero at the destabilization temperature T ~ 93 GeV. 
So all the phase transition happens in the 3 GeV interval shown in Fig. 1. Using that 
S^{T C )/T C ~ 140 — 145 one can easily compute the value of the temperature T c For 
the case considered in Fig. 1 we obtain T c ~ 95.2 GeV. The corresponding profiles 
for the Higgs bubbles h\{r) and h 2 (r) are plotted in Figs. 2 and 3 [where p 2 (r) = 
h\{r) + hl{r) and tan/3(r) = h 2 {r)/h x (r), A(3(r) = (3(r) - (3(T), (3(T) being defined by 
P(T) = v 2 {T)/v 1 {T)), for T = T c (solid lines), T = T c + 0.4 GeV (long-dashed lines) 
and T = T c — 0.4 GeV (short-dashed lines). 

The variation of the bubble parameters with respect to rriA is displayed in Figs. 4, 5 
and 6. In Fig. 4 the profile p{r)/v{T c ) is plotted for m A = 100, 200, 300 and 400 GeV, 
thick-solid, long-dashed, short-dashed and thin-solid curves, respectively. In all cases 
the bubbles have thick walls, the wall thickness being 

L u ~ (20 - 30) /T c 

as can be seen from Fig. 4. We can also see that wall profiles are almost indistinguisable 
for rriA ~ 200 GeV. In Fig. 5 we plot the parameter A/3 = A/3(oo) — A/3(0) as a function 
of tua for the whole effective potential, including the two-loop thermal corrections (solid 
line) and excluding them -i.e. including only the one-loop thermal corrections- (dashed 
line) . We see there is an enhancement due to the two-loop corrections which goes from 
~ 5 for rriA ~ 100 GeV, to ~ 2 for tua ~ 400 GeV. Since the total amount of pro- 
duced baryon asymmetry is proportional to Af3, the latter enhancement translates into 
a baryogenesis enhancement, which was disregarded in previous analyses |p8|| . Further- 
more, the enhancement produced in the strength of the phase transition, v(T)/T, by 
the two- loop thermal corrections [p5||-||28|| also affects the amount of baryon asymmetry 



produced, since the latter is proportional to [[UJ the integral 



DC 







dT T2 dr ' {6A) 



3 The actual temperature T c at which the transition happens is readily computed by comparing 
the probability of bubble nucleation per unit time and unit volume, ~ T A exp{— Ss(T)/T}, with the 
expansion rate of the universe at the corresponding temperature, r _1 = T 2 /AIp\, with = 
4tt^/tt [gsiT) + 7/8 gp(T)] /45. By imposing the condition that the probability for a single bubble to 
be nucleated within one horizon volume is 0(1) one can compute the temperature T c , i.e. from 

4 
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and therefore the two-loop enhancement is further strengthened. In Fig. 6 we plot the 
integral ( |3.1|) as a function of tua, when the two-loop thermal corrections are included 
(solid line) and excluded (dashed line). In all cases the total enhancement is greater 
than one order of magnitude. 



4 Squark bubbles 



In the previous analysis we have assumed that the bubble is built by the neutral field 
Higgses. This ansatz, used in the Euler-Lagrange equations, is obviously consistent. 
However, if we want to interpret this solution as the one controlling the tunneling 
process, we have to analyze the variation of the Euclidean action under generic quadratic 
perturbations including the rest of the fields. In particular, we have to make sure that 
there is just one negative mode: the breathing mode. 

The main candidate for such a new negative mode is provided by squark fields, in 
particular stop fields. If we include these field in our analysis, then the finite temper- 
ature effective potential will be now a function of four fields, V e ff(hi,h 2 ,tL,t R ). In 
general, there will be directions involving the new fields ti,t R where the potential de- 
creases []. Since the kinetic term is always positive, negative modes imply necessarily 
negative values for the variation of the potential energy. In the quadratic approxima- 
tion, this variation is given by 



<5Vr=47r r 2 dr 



t* L (r) t* R (r) 



M- 



tL(r) 
tn(r) 



(4.1) 



where 



M\(h 1 ,h 2 ) 



m 2 Q + U lL + h 2 t h 2 2 (r) 



A t h t h 2 {r) - fi hth\(r) 



A h t h 2 (r) - fj, h t hx{r) 



m, 



+ U 7R + h 2 t hl' 



(4.2) 



and n^- stand for the leading T— dependent self-energy contributions to the thermal 
masses. They are proportional to T 2 |19|, |34| . The eigenvalues of M. ~ (hi, h 2 ) give us the 
masses of the stop squarks in the Higgs background. We have to evaluate this masses 
along the path in the (hi, h 2 ) plane described by the bubble when r changes from zero 
to infinity. Since the variation of tan/5 along this path is less than O(10 -2 ), we will 
assume that A t is constant. Notice that these masses must be positive at the false and 
at the true vacuum because these are a local and a global minimum, respectively. 



4 This is certainly the case for some values of the parameters where new minima do appear. 
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These conditions translate into: 



ml + Ii TL > 



m 



+ n ?B > (4.3) 



( m Q + U t L +™ 2 t)(™ 2 U + ^l R + m t) > A t m t 

On the other hand, a necessary condition for having negative masses along the path 
followed by the bubble is: 

A 2 t >m 2 Q + U TL +m 2 u + U TR (4.4) 
However, this condition is forbidden by the requirement of not wipping off, after the 



phase transition, any previously generated baryon asymmetry ||28|| , and is never satisfied 
in our choice of parameter space. We then conclude that the bubble solution obtained 
in section 2 is not disturbed by non trivial configurations for the fields and/or 
and our conclusions in section 3 are rather robust. Had we chosen a value of the mixing 
A 2 satisfying condition ( |4.4j ) we could have in the bubble wall a field configuration with 
non-trivial values of £r and/or tt explicitly violating baryon number. 



5 Conclusions 

In this paper we have computed the tunneling processes from the symmetric phase to 
the true minimum in the first order phase transition of the Higgs fields in the MSSM. We 
have obtained the corresponding Higgs profiles along the bubbles. We paid particular 
attention to the so-called light stop scenario, where the phase transition is strong enough 
not to wipe off any previously generated baryon asymmetry. Baryogenesis in the MSSM 
was previously proved to be controlled by some bubble parameters, and in particular 
by the integral (|3.1|) , which were estimated by some energetic considerations based on 
the one-loop effective potential. We have shown, using our numerical calculations based 
on the two-loop effective potential, that there are important enhancement effects, with 
respect to those estimates, that can be around one order of magnitude. Finally we have 
proved that, for the considered cases, our solution is not disturbed by any non-trivial 
configurations of the stop fields. 
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Figure 1: The Euclidean action as a function of the temperature for m Q = 1 TeV, 
tan/5 = 2.5, A t = 0,m T = 150 GeV, and m A = 200 GeV. 
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Figure 2: The Higgs profile p(r)/v(T) for T = T c + 0.4 GeV (long-dashed curve) and 
T = T c — 0.4 GeV (short-dashed curve), and values of supersymmetric parameters as 
in Fig. 1. 
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Figure 3: The same as in Fig. 2 but for the Higgs profile A/5(r). 
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Figure 4: Higgs profile p(r)/v(T c ) for itlq = 1 TeV, tan/3 = 2.5, A t = 0, mj = 150 
GeV, and vtla = 100 GeV (thick-solid curve), 200 GeV (long-dashed curve), 300 GeV 
(short-dashed curve) and 400 GeV (thin-solid curve). 
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Figure 5: The parameter A(3 in the two- loop (solid curve) and one-loop (dashed curve) 
approximations, for mg = 1 TeV, tan/5 = 2.5, A t = and mj= 150 GeV, as a function 
of rriA- 



16 




Q I I I I I I I I I I I I 

100 200 300 400 

m A (GeV) 

Figure 6: The same as in Fig. 5 but for the integral U (3.1), £=l-loop (dashed curve), 
£=2-loop (solid curve). 
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